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Abstract 

Magnetic field and electrostatically tuned superconductor-insulator (SI) transitions of ultrathin metal films with levels of 
disorder that place them near the disorder-tuned SI transition appear to be direct, continuous quantum phase transitions. 
When films with lower levels of disorder are subjected to a perpendicular magnetic field, instead of a direct transition, a mixed 
superconductor-nonsuperconductor regime emerges at the lowest temperatures. The zero temperature limit of the resistance is 
either insulating or superconducting, depending upon the value of the field, suggesting that the behavior in this limit is governed 
by percolation physics. At high fields and low temperatures, in the nominally insulating regime, the resistance rather than 
the conductance is found to be a logarithmic function of temperature corresponding to predicitons for the weak localization of 
bosons. 



Highly-disordered, homogeneous, quench-condensed, 
ultrathin films of metals can be tuned between supercon- 
ducting and insulating behavior by magnetic field, thick- 
ness or carrier density [TJ [2]. These transitions are of 
interest because they are among the simplest of quan- 
tum phase transitions, and the field and carrier density 
tuned transitions appear to belong to the (2+l)D XY 
universlity class. Recently there has been increased in- 
terest in highly-disordered, amorphous compounds such 
as R12O3 and polycrystalline TiN, ranging in thickness 
from tens to hundreds of Angstroms. In addition to un- 
dergoing superconductor-insulator (SI) transitions tuned 
by perpendicular field, these systems exhibit magnetore- 
sistance, R(B), peaks in the insulating regime|3H7|. Al- 
though this effect was first observed and interpreted 
as evidence of a Bose insulating regime by Paalanen, 
Hebard, and Ruel [5] almost two decades ago, recent ex- 
periments have displayed enhancements of resistance by 
many more orders of magnitude j9j [10]. Furthermore, for 
less disordered films, an intermediate metallic phase has 
been reported to occur between the superconducting and 
insulating regimes [TT]. What is surprising is that in work 
on the much thinner, atomically disordered (amorphous), 
quench-condensed films, peaks in R(B) in the insulating 
regime of have only been observed in films patterned with 
a nanohoneycomb array of holes[12 , and there appears 
to not be an instrinsic intermediate metallic regime. In 
most of the measurements, normal resistances were close 
to the critical values associated with thickness-tuned SI 
transition. The present study, which was motivated by a 
search for magnetoresistance peaks in quench-condensed 
films, extends the regime of parameter space to normal 
resistances below criticality, yet not so far below that the 
field-tuned transition is to a metallic state. 

There have been two striking findings. The first is 
the behavior in the high-field limit, where the resistance, 
R(T), rather than the conductance, G(T), was found to 
vary as InT. This is suggestive of the weak boson localiza- 



tion phenomenon predicted by Das and Doniach [13], and 
reported for under-doped cuprates whose superconduc- 
tivity was completely quenched by magnetic fields |14j . 
In addition there was no peak in R(B) in fields up to 10 
T. A second finding is that at high temperatures (around 
2K ~ 4K) the data could be collapsed using a finite size 
scaling analysis. With decreasing temperature scaling 
broke down because curves of R(T) became nonmono- 
tonic functions of temperature in the manner suggestive 
of thickness tuned SI transitions of granular films |15j . 
This change in the physics with decreasing temperature 
is reminiscent of the approach to a quantum critical point 
of some stongly correlated electron systems in a magnetic 
field, in which new physics turns on before the quantum 
critical point is reached [15] • Mixed phases in the case of 
the field-driven SI transition have been discussed in the 
theoretical literature [171 \TE\ . 

The present investigations were carried out using films 
grown on (100) SrTiC>3 (STO) single-crystal substrates. 
Platinum electrodes, 100A thick, were deposited ex situ 
onto the substrate's epi-polished front surfaces to form 
a configuration suitable for both resistance and Hall re- 
sistance measurements. The substrate was then placed 
in a Kelvinox-400 dilution refrigerator/UHV deposition 
apparatus [H]. A 10 A thick under-layer of a-Sb and suc- 
cessive layers of a-Bi were thermally deposited in situ un- 
der ultra-high vacuum conditions through shadow masks 
onto the substrate's front surface. The substrate was 
held at about 7 K during the depositions. Films grown 
in this manner are believed to be homogeneously disor- 
dered on a microscopic, rather than on a mesoscopic scale 
[2U] . The sample measurement lines were heavily filtered 
so as to minimize the electromagnetic noise environment 
of the film. To avoid complications arising from this fil- 
tering, measurements were made using DC, rather than 
AC, methods. The ion and turbo pumps connected to 
the growth chamber were electrically isolated from the 
cryostat using vacuum nipples with ceramic spacers. 
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Figure 1: a-Bi (with a 10A thick a-Sb underlayer) film in 
different perpendicular magnetic fields. The values of field 
from top to bottom are 10, 5, 2.5, 2, 1.8, 1.6, 1.4, 1.3, 1.2, 
1, 0.8, 0.6, Tesla. Inset: Sheet resistance vs. temperature 
for the thickness-tuned superconductor-insulator transition of 
granular (without an underlayer) Bi films adopted from Ref. 
15. The thicknesses are 21.11 (top), 21.41, 21.75, 21.83, 21.94, 
22.04, 22.10, 22.14, 22.17, 22.23, 22.38, 22.63, 22.95, 23.37, 
23.66, 24.52, 25.40, 26.07, and 27.56 A(bottom). Notice that 
these two plots involve films with different morphologies and 
different tuning methods. 
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Figure 1 shows a series of curves of R(T) as a function 
of magnetic field charting the transition from supercon- 
ductor to insulator. 

The significant feature is the non-monotonic behavior 
of R(T) at values of magnetic field that are close to those 
separating superconducting and nonsuperconducting be- 
havior. This strongly suggests that there is a range of 
magnetic field over which the system breaks into a mix- 
ture of superconducting and nonsuperconducting regions 
at nonzero temperatures despite the homogeneity of the 
film. Behavior of this sort is usually found for granular 
films. To make the point, the inset of Fig. 1 shows the 
evolution from insulator to superconductor as a function 
of thickness for a granular film, which also shows a non- 
monotonic variation of R(T) |15| . but with much larger 
scale in both resistance and temperature. 

This 11. 2 A thick film is homogeneous even though it 
exhibits non-monotonic behavior at nonzero temperature 
within a range of fields over which it is undergoing a 
transition from superconducting to insulating behavior 
in the low temerature limit. Hall effect measurements 
support this. They reveal the areal carrier concentration 
to 1.5 x 10 16 cm~ 2 , which is close the value expected for a 
metal. From the carrier concentration, one can obtain the 
Fermi wavevector kp =2.5 x 10~ 8 cm -1 . Combining this 
with result for the sheet resistance, which gives kpi ^ 6, 



Figure 2: G(T) in zero magnetic field, (b) Scaling analysis 
for the perpendicular magnetic field tuned transistion between 
2K < T < AK and 0.6T < B < 5T. Resistance, normalized 
to the critical resistance (338457), is plotted as a function of 
\B - B c \/T 1/vz with vz=0A and B C =1.58T. 

we find the electronic mean free path t to be around 2.4A. 
This indicates that the disorder is on an atomic rather 
than a mesoscopic scale. The carrier concentration would 
not be expected to be as high in a granular film. 

We now turn to the details of the transport. The con- 
ductance G can be fit by InT in the normal state in zero 
field and in low magnetic fields. The coefficient of InT is 
5.23 x 10~ 6 f2 _1 D which agrees with the theory of weak 
localization including electron-electron interactions, for 
the case of strong spin-orbit coupling |2lj . The coefficient 
of the logarithm also agrees with that found in previous 
work [2]. A sample fit, with the range of applicability 
clearly delineated is shown in Fig. 2a. 

In the vicinity of the crossover magnetic field, the re- 
sistance at relatively high temperatures could be loosely 
fit with the form R ~ Rexp(To/T), with a negative value 
of T for films that smoothly become superconducting 
and a positive value for those which ultimately become 
insulators. This suggested that R(T, B) might be col- 
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lapsed using the finite-size scaling relation suggested by 
Fisher [22] • Figure 2b shows the result of this analysis for 
2K < T < AK and 0.6T < B < 5T employing the form 
R/R c ~ F (\B - B^/T 1 /"*) with the critical exponent 
product vz — 0.4. The scaling fails when applied to data 
outside of these ranges. 

Well into the insulating regime, where one might ex- 
pect a return to G ^ InT, but one finds instead, that R 
rather than G is better described by InT. A comparison 
of these two forms for data obtained at a field of 10T is 
presented in Fig 3. This relationship was also found at 
lower values of magnetic field. We compared the values of 
chi square x 2 = J2 [(Ri ~ Rfit{Ti)) / <Ji\ for the two func- 
tional forms (R ^ InT and G ^ InT), where the Ri are 
the data points, Rf it {Ti) is the fitted function, and the 
Oi are the errors of the data points. We include data over 
the range from 100 mK to 1 K, and in magnetic fields of 
2.5, 3, 4 and 5 T. For R - InT the values were 6104, 7063, 
2932, and 1291 respectively, whereas for G ^ InT they 
were 37276, 40298, 19957, and 8216, respectively. This 
makes the case for R ^ InT being the better description 
of the data. 

Futhermore, if one forces the data to be fit by the func- 
tional form G ^ InT, as shown in Fig 3(b), the coefficient 
of InT is nearly three times larger than that found in the 
high-temperature low-field regime. For example, in zero 
field in the temperature range of 5K ~ 10K , the coeffi- 
cient is 5.23 x lO^rr 1 ::, whereas it is 1.63 x lO^O^D 
at 10T over the temperature range from Q.IK ~ IK 
when forced. This indicates that the resistance of the 
film in the high field, low temperature regime increases 
more rapidly than that of a quantum corrected metal. 

Figure 4 provides a global representation of the regimes 
over which R(T) can be described by various forms as 
well as crossover and nonmonotonic regimes which can- 
not be described by simple functional forms. The data 
clearly suggest that the evolution from a superconducting 
ground state to a nonsuperconducting state is indirect, 
involving an intermediate regime at nonzero temperature 
in which there are superconducting islands in a nonsuper- 
conducting matrix[T7]|18|. 

In the case of quench-condensed, disordered supercon- 
ducting films with disorder levels closer to the disorder- 
or electric field-tuned transitions than the present films, 
in the high-field regime the behavior is clearly insulating 
although there is no magnetoresistance peak which has 
been taken as evidence that the insulator is a Cooper pair 
insulator. The latter behavior has only been observed 
in compound amorphous films or in quench condensed 
films patterned with nanohoneycomb hole arrays. The 
high magnetic field regimes of the films of the present 
work are not a quantum corrected 2D metals, because 
as demonstrated above, R(T) rather than G(T) is better 
fit by InT over an extended range of temperature. The 
temperature dependence is the same as that predicted in 
the work of Das and Doniach which describes a Bo- 
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Figure 3: logT with the results of fit in a 10 T field. The red 
thin lines are the fitting lines. The resistances are fit better 
by InT over a larger temperature range (up to 6K) than the 
conductances. 



son weak localization model. Although the coefficient of 
the InT term is an order of magnitude smaller than that 
which they predict, the model may be relevant. Similar 
behavior has been observed in underdoped cuprate su- 
perconductors when superconductivity is quenched by a 
magnetic field [14] . Moreover, the coefficient of InT in 
the present work is close to that found in that work. If 
one converts the sheet resistance into resistivity, the co- 
efficient of InT in a field of 2T is —6.9 x 10~ 5 Q-cm, while 
for the case of x=0.13 of La2-xSr x Cu04, it was found to 
be —7.2 x 10~ 5 f2-cm when current is in the ab plane. 

The fact that at higher temperatures the data can be 
scaled, suggests the existence of a quantum critical point. 
However, a direct quantum phase transition never occurs, 
because the two-phase regime develops as one lowers the 
temperature. In the limit of zero temperature the film 
is either insulating or superconducting so that it is likely 
that percolation physics plays a role in the SI transi- 
tion. The value of the critical exponent product vz = 0.4 
found in the regime in which the data scales is not as- 
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Figure 4: R ~ (?iT and G ~ InT regimes correspond to 0.1% 
difference between data points and the results of the fits. The 
boundaries of the scaling regime were set using the same cri- 
terion. The indicated value of B c corresponds to the critical 
field of the scaling analysis. 



sociated with any particular model relevant to ultrathin 
films. However this value has been found in a Monte 
Carlo simulation of (2+l)D XY Josephson junction ar- 
rays without disorder and with frustration / = 1/4 |23j . 
Furthermore, this simulation exhibits a first order phase 
transition when the frustration drops to f — 1/5. A two- 
phase regime such as inferred from our data could also 
imply a first order quantum phase transition. Model- 
ing a disordered film with a junction array would appear 
to require a Monte Carlo simulation including disorder, 
which is not included in the simulation of Lee and Cha. 
On the other hand, Kim and Stroud performed a simula- 
tion including disorder but without magnetic field [S3] . A 
simulation with disorder and frustration would be needed 
to model films such as those reported in this work. 

The present results deepen the mystery associated with 
the superconductor-insulator transition. Although vari- 
ous models may account for the behavior of specific sys- 
tems, there is no model that explains their different be- 
haviors. For example, it is not clear as to why some 
transitions are direct, and with small changes in the level 
of disorder an intermediate two-phase regime emerges. 
Also the reasons for the differences between the high field 
regimes of different materials are not known. 

This work was supported in part by the National Sci- 
ence Foundation under grants NSF/DMR-0455121 and 
NSF/DMR-0854752. 
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